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We establish a general thermodynamic scheme for cosmic fluids with internal self-interactions and 
discuss equilibrium and non-equilibrium aspects of such systems in connection with (generalized) 
symmetry properties of the cosmological dynamics. As an example we construct an exactly solvable 
gas dynamical model of a "deflationary" transition from an initial de Sitter phase to a subsequent 
Friedmann-Lemaitre-Robertson- Walker (FLRW) period. We demonstrate that this dynamics rep- 
resents a manifestation of a conformal symmetry of an "optical" metric, characterized by a specific 
effective refraction index of the cosmic medium. 



I. INTRODUCTION 

A comprehensive picture of the physics of the early 
universe crucially relies on an adequate understanding 
of thermodynamic aspects of the cosmological evolution 
jl],^| . Of particular interest are periods during which the 
cosmic substratum was out of equilibrium. It is gen- 
erally believed that essential properties of the presently 
observed universe may be traced back to specific non- 
equilibrium processes in early cosmology. A detailed and 
transparent description of out-of equilibrium phases is 
therefore a basic challenge for any cosmological model. 
On the other hand, the natural starting point for thermo- 
dynamical studies is a characterization of possible equi- 
librium states of the system under consideration. It is 
well known that equilibrium states in the expanding uni- 
verse are only possible under certain conditions. For a 
simple gas model e.g., the cosmic matter has to obey ex- 
actly an equation of state p = p/3 {p is the pressure and 
p is the energy density) |||| . This equilibrium condition 
is equivalent to a symmetry requirement for the cosmo- 
logical dynamics: The quantity u a /T, where u a is the 
four velocity of the medium and T is its temperature, 
has to be a conformal Killing vector (CKV). For devia- 
tions from p = p/3 the conformal symmetry is violated 
and no equilibrium is possible. As long as the deviations 
from equilibrium are small, there exist well-established 
approaches and techniques to investigate the correspond- 
ing dynamics (see, e.g. f|-|(|). For larger deviations the 
situation is less clear and new concepts are required. 

In this paper we focus on a specific class of non- 
equilibrium states which are characterized by a simple 
modification of the standard equilibrium and symmetry 
conditions for relativistic media. We demonstrate how 



the standard equilibrium of relativistic gas dynamics may 
be generalized such that ultrarelativistic matter is no 
longer a singular case but a well-defined limit of a more 
general equilibrium configuration. This generalization is 
achieved by giving up the simple dilute gas approxima- 
tion and admitting additional internal interactions in- 
side the many-particle system. As long as these inter- 
actions are particle number preserving the correspond- 
ing states are real equilibrium states. The step away 
from equilibrium is then performed by including internal 
interactions which do not preserve the particle number. 
This procedure, although formally rather straightforward 
within the presented formalism, is physically essentially 
non-trivial since it introduces states with (not necessar- 
ily small) entropy production. It is a remarkable fea- 
ture of this kind of non-equilibrium that, irrespective of 
the amount of entropy production, it is characterized by 
an equilibrium distribution ("generalized" equilibrium) 
of the microscopic particles and is accessible by an an- 
alytic treatment. Furthermore, there exists a threshold 
for the entropy production in homogeneous and isotropic 
models beyond which such type of states are only possible 
in a universe with accelerated expansion. This represents 
a surprizing link between (generalized) equilibrium prop- 
erties of cosmic matter and the expansion behaviour of 
the universe. 

The present paper generalizes and extends previous in- 
vestigations ]7|-[lCt| on specific aspects of generalized equi- 
librium configurations in a cosmological context. Our 
first purpose is to provide a general characterization of 
all those (equilibrium and non-equilibrium) states of a 
gaseous cosmological fluid which on the microscopic level 
are governed by an equilibrium distribution function. 
The second purpose is to establish on this basis an exactly 



1 



solvable homogeneous and isotropic gas dynamical model 
according to which the universe starts in a phase of ex- 
ponential expansion and subsequently smoothly evolves 
towards a standard radiation dominated FLRW phase. 
This "deflationary" jllj transition is shown to be the 
manifestation of a specific non-equilibrium configuration 
of the cosmic medium. We determine those forces on the 
microscopic constituents of the cosmic matter which un- 
derly this dynamics. Moreover, we demonstrate that a 
non-equilibrium transition of this kind represents a con- 
formal symmetry of an "optical" metric |12-14| with an 
initially large value of the "refraction" index of the uni- 
verse which continously decreases to unity as the FLRW 
phase is approached. 

The paper is organized as follows. In Sec. II we es- 
tablish the general thermo-hydrodynamical framework 
for a relativistic fluid with particle production and ef- 
fective viscous pressure. This scheme essentially relies 
on a consistency requirement for the temperature evo- 
lution of the medium which relates symmetry aspects 
with those of thermodynamics. A corresponding micro- 
scopic realization is given in section III, where we study 
the kinetic theory of a gas of self-interacting particles. 
This encompasses a determination of the effective self- 
interacting force by generalized equilibrium requirements 
and a derivation of corresponding transport equations, 
followed by an investigation of the microscopic particle 
motion in such a system. In Sec. IV we demontrate how 
a specific self-interaction in a homogeneous and isotropic 
cosmic medium implies a deflationary (Tl| l transition from 
an initial de Sitter phase to a subsequent FLRW period. 
We determine an optical metric with respect to wich the 
quantity u a /T is a CKV for the entire transition pro- 
cess. An equivalent two-component description is given 
in Sec. V, where the deflationary dynamics is interpreted 
in terms of a decaying vacuum. Section VI sums up the 
main features and conclusions of the paper. Units have 
been chosen so that c= ks =Ti= 1. 



II. THERMO-HYDRODYNAMIC APPROACH 

We assume the cosmic substratum to be modeled by a 
fluid with the energy-momentum tensor 
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where p is the energy density measured by an observer 
comoving with the fluid four- velocity it 1 , normalized by 
u a u a = — 1, and hik = 9ik + UiUk- The total pressure P 
is the sum 
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of an equilibrium part p > and an effective non- 
equilibrium part LI < which is connected with entropy 
production. Within the Eckart frame the particle num- 
ber flow vector N l is given by 



where n is the particle number density. While we require 
local energy-momentum conservation T^^y k = in ac- 
cordance with the integrability conditions of Einstein's 
equations, we will not assume the fluid particle number 
to be conserved. Instead, we admit a source term V in 
the corresponding balance equation which describes the 
rate of change of the fluid particles: 



AT* = fi + On 



nT 



(4) 



where = u\ is the fluid expansion and n = n^u % . For 
r > we have particle creation, for T < particles are 
annihilated. The particle number may change either due 
to internal reactions within the medium [ |15|flq ] or due 
to particle production in strongly varying gravitational 
fields [p7|-p3| . The production rate T is an input quantity 
in a phcnomcnological description which in principle has 
to be calculated from the underlying microphysics. In the 
general part of this paper T will be a free parameter. In 
the subsequent application a specific dependence on the 
Hubble rate of an homogeneous and isotropc universe is 
considered. 

Energy-momentum conservation is equivalent to 



and 



p + e(p + p + Tl) = 



(p + P + n) ii a + v a (p + n) = o 



(5) 
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where ii a = u a . c u c is the four-acceleration, and V Q p = 
h h a pfi is the covariantly defined spatial derivative. With 
the help of the Gibbs equation (see, e.g., ||) 



Tds = d— + pd— . 

n n 



(J) 



where s is the entropy per particle and T is the equilib- 
rium temperature, we obtain 

Ti 

nTs = p-(p+p)- . (8) 

n 

Using here the balances (||) and @ yields 

nTs = -en- (p+p)T . (9) 

From the Gibbs-Duhem relation (see, e.g., ||) 

dp=(p + p)^+nTd^) , (10) 
where p is the chemical potential, it follows that 
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p p + pT 



(11) 



.77 nT nT T 
We assume equations of state in the general form 

p = p(n,T) , p = p(n,T) , (12) 
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i.e., particle number density and temperature are taken 
as the independent thermodynamical variables. Differen- 
tiating the latter relation and using the balances (0) and 
(|J) provides us with an evolution law for the tempera- 
ture: 



T [ 1 dp dp/dT ' 



(13) 



where the abbreviations 

dp _ (dp/dT) n 



dp 
df 



dp 
df 



dp ~ (dp/dT) n ' 
have been used, as well as the general relation 
p + p 



dp 
dn 



T dp 
ndf ' 



which follows from the fact that the entropy is a state 
function, i.e., d 2 s/dndT = d 2 s/dTdn. With the help of 
Eqs. (|l^), (|]), (||) and ( |l3| ) p may generally be written 
as 



c 2 s P 
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where 



dp\ n dp T (dp/dT) 2 

dp) ad p + pdn p + p dp/dT 
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is the square of the adiabatic sound velocity c s [see, e.g. 
M], Using Eq. (El) in Eq. (Til) yields 
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To obtain an expression for the spatial variation of p,/T 
we combine the Gibbs-Duhem equation ( |To| ) with the mo- 
mentum balance (^). The result is 
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Together with the particle number balance (|J) the en- 
tropy flow vector S a = nsu a gives rise to the following 
expression for the entropy production density: 



S° a = nsT- 



S a = nsu a = 
It is useful to rewrite the latter as 
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where £^ denotes the Lie derivative with respect to the 
temperature vector £ a = u a /T . The equivalence between 
s in Eq. (||) and the right-hand side of the last relation 
may be checked by direct calculation. The purpose of 
this rewriting is to relate the entropy production density 
to the Lie derivative of the metric which, in turn, may 
be used to characterize symmetries of the spacetime. 

As any symmetric tensor, the Lie derivative may be 
split into contributions parallel and perpendicular to the 
four- velocity: 

£^9ik = 2AuiU k + B x u k + B k u t + 2<ph lk + b ik , (20) 
where 
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and 
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the quantity o a b being the shear. 

It is crucial for the following considerations that the 
temperature behavior of the fluid is governed by the gen- 
eral thermodynamics of the fluid according to the law 
( |l3| ) and, on the other hand, by the projection ( pl| ) of 
the Lie derivative. Consistency between the expressions 
(H) and (|l|) requires 
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n s 



dp 

dp~ + 2 Tdp/dT 
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All relations considered so far, including the last one, be- 
come especially simple for s = 0. Under this condition 
the entropy production is entirely determined by the par- 
ticle production rate, i.e., S a a = nsT. For a conserved 
particle number the entropy production vanishes. Only 
due to the enlargement of its phase space the system in- 
creases its entropy content, i.e., each newly created parti- 
cle contributes with its (equilibrium) entropy to the total 
entropy of the system. There are, however, no dissipa- 
tive processes connected with a change in the entropy per 
particle. Except for a change in the particle number (of 
whatsoever reason) the system behaves as a perfect fluid. 
According to Eq. (^|) the condition s — establishes a 
simple relation between the non-equilibrium pressure II 
and the particle creation rate F [|25|]2£ 
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s = => U = -(p + p)-. (26) 

With such type of viscous pressure the energy balance 
(||) is equivalent to 

p + e(p + p) = (p + p)T . (27) 

Consequently, the system may be understood as a per- 
fect fluid with source terms. The microscopic counter- 
part (see below) of this statement is the fact that corre- 
sponding states for a gas are described by an equilibrium 
distribution function. 

With relation (|2^) and for s — 0, the general structure 
of the Lie derivative becomes 
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We emphasize that the consistency requirement ( J25| ) 
(with s = 0) on which the expression ( |28| ) is based only 
affects the projection u l u k £^gik but not the parts in- 
volving ui + ViT/T and er^. The introduction of the 
parameter A is convenient since it allows us to handle 
the cases of preserved and non-preserved particle num- 
bers simultaneously. Particle number conservation T = 
is equivalent to A = 1. Particle production in Eq. J2S| ) 
and all further relations involving A is then straightfor- 
wardly taken into account by substituting A = 1 by the 
quantity (p9|). As already mentioned, this formally sim- 
ple replacement implies the transition from equilibrium 
configurations with vanishing entropy production to non- 
equilibrium situations with S? a > 0. In the limit Oik = 0, 
iii + V. L T/T = 0, A = 1, and dp/ dp = 1/3 Eq. (§|) 
reduces to the CKV equation for u l /T, i.e., to the famil- 
iar global equilibrium. In the following we briefly list up 
other previously investigated special cases of Eq. ( p8| ) : 

(i) A simple extension of the CKV condition is ob- 
tained by keeping — 0, iii + ^iT/T — 0, A = 1, but 
allowing dp/ dp to be different from 1/3. The CKV prop- 
erty is modified in such a case, but there is no entropy 
production, i.e., the corresponding states are equilibrium 
states (l(J. This represents an extension of the global 
equilibrium concept to equations of state different from 
p = p/3. Ultrarelativistic matter is then the limit in 
which the condition ( p8| ) (with A = 1) reduces to the 
CKV condition for £ a . 

(ii) A special case which implies a transition to non- 
equilibrium states corresponds to the vanishing of the 
parenthesis in front of u^k within the brackets in the 
first term on the right-hand side of Eq. (EM) : 
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If, moreover, in + ViT/T = and cr^ = are valid, we 
have £^gik = 24>9ik, i-e., u % /T is a CKV. For A = 1, cor- 
responding to r = 0, i.e. particle number conservation, 
the condition (|3^) necessarily implies the ultrarelativistic 
equation of state p = p/3. For r / 0, however, equiva- 
lent to A / 1, this condition is compatible with different 
equations of state as well [0J|]. Admitting r ^ corre- 
sponds to a non-equilibrium generalization of the global 
equilibrium. Again, ultrarelativistic matter is then no 
longer a singular case but just the limit for r — > 0. 

(iii) Another interesting case is T = equivalent to 
A = 0, which corresponds to a stationary temperature 
behavior f = 0. For u t + V t T/T = and a lk = the 
CKV condition £^gik — 2(f>gik is replaced by £^gik — 
2<t>fkk H ■ In Sec. IV the latter relation will be used to 
characterize a de Sitter phase in the early universe. 

It is obvious that there exist further generalizations of 
the global equilibrium which in general both modify the 
conformal symmetry and imply particle production, i.e. 
A / 1. The CKV condition £^gik = 2(pgik is known 
to follow from Boltzmann's equation as global equilib- 
rium condition for the one-particle distribution function 
of a classical ultrarelativistic gas. For the special cases 
mentioned above it was shown that the corresponding 
modifications of the standard conformal symmetry with 
p = p/3 naturally appear as "generalized" equilibrium 
conditions for the one-particle distribution function for 
particles which self-consistently move in additional force 
fields [p|-^0[. In the follo wing section we extend these 
results to the general case (|2g|). 



III. KINETIC THEORY IN A FORCE FIELD 

A. Generalized equlibrium conditions 

The particles of a relativistic gas are assumed to move 
under the influence of a four-force F l in between elastic 
binary collisions, described by Boltzmann's collision inte- 
gral J [/, /] . The equations of motion of the gas particles 
are 



da; 1 
d7 
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Tip 1 
fry 



= mF l 



(31) 



where 7 is a parameter along the particle worldline which 
for massive particles may be related to the proper time 
t by 7 = r/m. Since the particle four-momenta are nor- 
malized according to p l pi = — to 2 , the force F % has to 
satisfy the relation piF 1 = 0. The corresponding equa- 
tion for the invariant one-particle distribution function 
/ = / (x (r) ,p (t)) may be written as 



dx % 



dp k 



dp 1 



Both the collision integral J and the force F l describe 
interactions within the many-particle system. While J 
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accounts for elastic binary collisions, F l is to model dif- 
ferent kinds of interactions in a simple manner. The 
force F l is not a fundamental force. Strictly speaking, 
F l should be calculated from the microscopic particle dy- 
namics and, consequently, depend on the entire set of par- 
ticle coordinates and momenta characterizing the system 
of gas particles. It is well known that except for simple 
cases a relativistic statistics for interacting many-particle 
systems is not available. In such a situation it seems 
natural to resort to a semi-phenomenological or effective 
theory. We assume that all the microscopic interactions, 
however involved their detailed structure may be, can 
be mapped onto an effective one-particle quantity, which 
may be used for a subsequent self-consistent treatment. 
The situation here is more complicated than for exam- 
ple for charged particles in an electromagnetic field were 
the microscopic starting point for a self-consistent anal- 
ysis is a given interaction, the Lorentz force. Here, the 
analytic structure of the force is not known in advance. 
We restrict ourselves to the class of forces which admit 
solutions of Eq. j3^ ) that are of the type of Jiittner's 
distribution function 



f(x,p) = exp{a + p a p a } , 



(33) 



where a = a (x) and j3 a (x) is timelike. For / — > /° the 
collision integral vanishes, i.e., J [/°, /°] = 0. Substitut- 
ing /° into Eq. ([32]) we obtain 

p a a a + P(a-b)P a p b = -mfrF 1 . (34) 

Generally, the effective force F l is an arbitrary function 
of the particle momentum. However, focusing on forces 
which satisfy the equilibrium condition ( p4] ) implies con- 
siderable simplifications. This becomes obvious if we ex- 
pand the quantity (3iF l on the right-hand side of Eq. (|34|) 
in a power series in p a : 

PiF 1 = AF l (x) + faFl (x)p a + ft^ 6 (x)p a p b + .... 

+/3 i F^ n (x)p a p b ....p n + ... (35) 

Here we assume PiF^ also to contain all those trace parts 
of odd higher-order terms (third order and higher) which 



by the condition pip 1 — 



may be reduced to first 



order in p a . Analogously, fiiF l ah is supposed to contain 
all those trace parts of higher even order (fourth order 
and higher) terms which by the condition pip 1 = —m 2 
may be reduced to second order in p a . This is equiva- 
lent to the assumption that the corresponding traces in 
PiF^ b (x) p a p b ....p n vanish for the third and all higher- 
order terms on the right-hand side of the expansion (35). 
We emphasize that it is only the projection faF 1 of F' 1 
which is expanded in a power series, not the force it- 
self. For the latter there exists the additional require- 
ment p a F a — 0. 

The point is now that the left-hand side of Eq. j3^ ) 
contains terms linear and quadratic in the particle mo- 
mentum. Consequently, the right-hand side of this equa- 
tion is restricted to quadratic order as well. It follows 
that the most general structure of the force term is 



PiF i = &F* (x) p a + fcFi b (x) p a p b , 
where we have used the redefinition 
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b kl o9kt 



Ft 



(36) 



(37) 



Comparing now different powers in Eq. (|34j) separately, 
we obtain the (generalized) equilibrium conditions 



and 



0(k;l) = -^-Epgki 



-mPiFl 



(38) 



(39) 



The last equation relates the Lie derivative of the metric 
to the force exerted on the particle. With the identifi- 
cation Pi = Ui/T we may read off the coefficient UiF^ b 
which is necessary to recover condition (Eq): 



mu t Fa b = j 



9ab 



Ua 

' 2 



u b 



1 3A^ ) a„m, + n nh 
dp) 

V b Tl u b \ . V a T 

u„ H 

T 2 T 



(40) 



This condition fixes the part which is quadratic in the 
particle momenta. The expression (J4(i|) is the most gen- 
eral quadratic force coefficient in Eq. (pq) which is com- 
patible with an equilibrium distribution function ([33]) . 
With a = /i/T the linear part is determined by 



m 



tu,f: = - 



(f)"«. + v.(£) • (41) 



An example for a force linear in p a is the Lorentz force 
which is obtained for F ab — > ^F^ b m ^ where e is the charge 
and F^ b m ^ is the electromagnetic field strength tensor. 
For this case the equilibrium condition ( |38|) is well known 
in the literature Other cases have been discussed 

in a cosmological contex t . 

With Eqs. ([39]) and (php we have recovered relation 
( p8| ) as equilibrium condition for a gas of particles mov- 
ing in the field of a force (|36|). It is a remarkable feature 
of this force that it depends on the kinematic quantities 
expansion, shear, acceleration and on the temperature 
gradient which are variables characterizing the system on 
a macroscopic level. These macroscopic quantities deter- 
mine the motion of the individual microscopic particles 
which themselves are the constituents of the medium. In 
other words, our system is a self-interacting gas. The 
analogy to the well known self-consistent coupling of the 
particle momentum of a charged particle to the electro- 
magnetic field strength tensor is obvious. The coeffi- 
cient UiFl corresponding to the linear part of the force 
is determined by the spacetime behavior of the macro- 
scopic variable /z/T, which will be relevant for the energy- 
momentum balance to be discussed below. 
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The relationship between Eq. ( J28| ) and Eqs. ( |39| ) 
and (^o|) admits an alternative interpretation which ex- 
hibits features that are familiar from gauge field theories. 
Gauge field theories rely on the fact that local symme- 
try requirements (local gauge invariance) necessarily im- 
ply the existence of additional interaction fields (gauge 
fields). In the present context we impose the "symme- 
try" requirement (^8|) (below we clarify in which sense 
the modification of the conformal symmetry is again a 
symmetry). Within the presented gas dynamical frame- 
work this "symmetry" can only be realized if one intro- 
duces additional interactions, here described by an effec- 
tive force field F % . Consequently, in a sense, this force 
field may be regarded as the analogue of gauge fields. 



B. Transport equations 

The particle number flow 4-vector N l and the energy 
momentum tensor T tk are defined in a standard way (see, 
e.g., |,1) as 



J dPp'f&p) , T ik = J dPp l p k f(x,p) . (42) 



The integrals in the definitions (|4^) and in the following 
are integrals over the entire mass shell p % p% = —m 2 . The 
entropy flow vector S a is given by 



S a = -Jp a [/In/-/] dP , 



(43) 



where we have restricted ourselves to the case of classical 
Maxwell-Boltzmann particles. The entropy production 
density is a sum of two terms: 



F 



Here. 



o c = - J dPJ[f, /] In/ 



(44) 



(45) 



is the familiar contribution of Boltzmann's collision inte- 
gral, while 



rrj, = ,„ / dPF'J^hl/ 



(46) 



takes into account an entropy production due to the ac- 
tion of the force F l . Since Boltzmann's H theorem guar- 
antees o~ c > 0, we have 



S?i -<j f >0 



(47) 



The equality sign in the last relation, realized by / — > f°, 
characterizes the generalized equilibrium as a state with 
minimal entropy production. With / replaced by / in 
the definitions © and @, the quantities N a , T ab and 



S a may be split with respect to the unique 4-velocity u a 
according to 

N a = nu a , T ab = pu a u b +ph ab , S a = nsu a . (48) 

Here we have identified the general fluid quantities of the 
previous section with those following from the dynam- 
ics of a Maxwell-Boltzmann gas. Note that the energy- 
momentum tensor T lk in Eqs. (ff2|) and (|4^) does not 
coincide with the tensor Tf^^ introduced in Eq. (Q). 
The exact integral expressions for n, p, p are (see, e.g., 



and 



P 

n 



p 4irm 2 T 
*i(f) 



K 2 (^) exp [a] , 



Mr) 



■ZT = m K ^\-T 

Mr) 



(49) 



(50) 



The quantities K n are modified Bessel functions of the 
second kind H . The entropy per particle s is given by 



The balances 

N a „ 



P + P _ £ 
nT T 



mft / F*/°dP , 



(51) 



T ak k = -m/3i I p a F i f°dP 



may be written as 

N« a = -mfr {F*N a + F* ab T ab ) 

and 

-mft {F l b T 



rpak 
1 -k 



i-p<tt> | -pi j^j-akl^ 



(52) 



(53) 



(54) 



where M akl = J dPf°p a p k p l is the third moment of the 
equilibrium distribution function. For the entropy pro- 
duction density we find 

S? a = mfa J[a + (3 aP a ] F i f°dP = -aN« a - (3 a T ab . 

(55) 

Usin g th e expressions ([hj|) and (|4l]) and the decomposi- 
tion (|4"§|), the particle number balance becomes 



(56) 



(£)' 


+ nO 


" P 






nT 



nT nT dp 



From Eq. ( J16| ) for s = (which implies relation (|26|)) we 
obtain 



p + p ( dp . , 



T 



dp 



c A9 



(57) 
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Since p = nT the expression (|l5|) simplifies to 
o nT 



1+^ 



nT 9p 



Ae 



Consequently, the balance (|56|) reduces to 
N? a = (1 - A) nQ = nF , 



(58) 



(59) 



which is consistent with Eq. (|j). For A = 1 the particle 
number is conserved. To evaluate the energy-momentum 
balance (|54|) one needs the relevant third moments of the 
distribution function, which may be obtained, e.g., with 
the help of the relation 



U aP a f = 



df° 

8 ft) 



(60) 



where the derivative has to be taken for a — const. One 
finds 



£~r(, + ,), ! = - 3 r ( , + P >- P r( ? 



(61) 



where we have used the properties 
d (K n {z)\ _ K n+1 (z) 



dz V z n 



K n 



, K n+ i = 2n h K n _i 

z 



(62) 



of the functions K n j|] . The relevant moments are 

g kl M akl = g kl f dPf P a p k p l = -nmV , (63) 



u a u kUl M aa = -3T (p + p) - pT ( - 



and 



Uk\Ul + 



T 



u a M akl = CJ k lU a M° 



(64) 



(65) 



For the energy balance we have 
u a T ak = -p{£j -{p + p)Q 

+nTQA— (3— + (- 
With the help of the relation (?J 



\t) 

it follows that 



(66) 



(67) 



Uo T^ = -(l-A)e(p + p) 



(68) 



The right-hand side vanishes only for A = 1. If we in- 
troduce the quantity II according to relations ( p6| ) and 
( ^9| ) we recover Eq. (|j). The relevant moments for the 



momentum balance h ca T a \ are 



Uk [ill + 



ViT 



T 



V C T 
T 



T{p + P) , 
(69) 



as well as UkUih ca M aM 
Eqs. (H), ©, © and @ one obtains 



0, and a k ih ca M 



akl 



0. From 



hcaT* = (p + p)u c + V c p 



nTV 



c (t) 



V C T 
T 



(p + p) • (70) 



In order to be consistent with Eq. (||) we have to have 
V T x 



nTV, 



T 



(P + P) 



-v m n - Uu r: 



(71) 



with the same II as above, given by Eq. (|26|). Conse- 
quently, the quantity V (fJi/T) is determined by 



-V m (|)=-^A 



(1-A) 



r 



r 



(72) 



This together with Eqs. (41) and (57) determines the the 
most general linear part —muiF l a p a of the force, compat- 
ible with an equilibrium distribution function (33). 



C. Microscopic dynamics 

Once the effective one-particle force is known, we may 
study the motion of the gas particles explicitly. We ex- 
tend here earlier investigations for special forces |8]-|Tc|l 
and consider again the most general case. Contracting 
the equation of motion in Eq. (|3l]) with the macroscopic 
four- velocity results in 



D ( M ^) rn 1 i k 

—± '- = u t F l + —u l , k p l p k , 

dr m 



(73) 



where we have used that 

Du l ■ p n 

-J- = U ;n— ■ 

dr ' m 

With the well-known decomposition of the covariant 
derivative of the four- velocity P7U28] , 



-UiU n + (J in + LU in + 



3 in 



(74) 
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where u) a t, — h^hfu] c;c n, we find 



IV. THE DEFLATIONARY GAS UNIVERSE 



D (u lP l ) 



dr 



= u. t F 
1 



3m 



Qh lk p l p k 
1 



i k • i k 

— (JikPP UiU k pp 

m m 



(75) 



The general structure of UiF % is given by Eqs. ( |36| 
and (Elf), 



u t F l = - 



(r) Ua 



e 

■g- I 9kl 



1-3A$ 



r 



1 /. , V;T 

1 /. V fc T 
-2 U '( U * + — 



U fc U; 



pV 



T 

rn 

_e_ 

3m 
E 



m \ J 



P k i 



(76) 



where B = —Uip 1 is the particle energy. Using the ex- 
pression J76| ) in Eq. ( [75|) together with Eq. jl5| ) for s = 
yields 



dr m 

T 



E T 
TT 



7 °(t) 



T T 



Since 



dr m m 



we obtain 



f ~ T 



= 



(78) 



(79) 



This relation has to be valid for arbitrary p m . It follows 
that 

_ „ u + u a p a E — Lt 
a + P a P = ^ = t^t- = const , (80) 



T 



T 



i.e., the equilibrium distribution ( p3\ j is indeed main- 
tained without any restriction for the equation of state. 
This proves the consistency of our approach in the gen- 
eral case. 



The cosmic medium is now regarded as a gas with 
internal self-interactions which macroscopically manifest 
themselves as an effective viscous pressure. For the dis- 
cussion of the cosmological dynamics we restrict ourselves 
to a homogeneous and isotropic spatially flat universe. 
Einstein's equations then reduce to 



8ttG P = 3H 2 



H = -4nG{p + p + Tl) , (81) 



where H = 8/3 = a/ a is the Hubble rate and a is the 
scale factor of the Robertson- Walker metric. Together 
with relation (1261) one finds 



r 

3H 



H 



3 (1 + w) H 2 



A 



H 



3(1 + w) H 2 



(82) 



where w = pip. According to Eqs. (|26j), (gg), (|36j), Qtf, 
(&T|), (p8|), (|72|), the viscous pressure is a consequence of 
the action of a force F % on the microscopic constituents 
of the cosmic substratum with a projection 



mu,F l 



7 



1 + w 
E 2 



w dp 



H 

H E 
dp H 



1 + w dp H 2 



H 



(83) 



For the previously discussed special case in which u l /T is 
a CKV with w — — < 1 and A = 1/2, Eq. @) specifies 

to a 



(77) 4% 



H 

H 2 



4 

H cx — 
3t 



oat 



4/3 



(84) 



We obtain accelerated expansion (a > 0) as a result of 
the back reaction of the production of massive particles 
on the Hubble parameter This dynamics is realized 
by a microscopic force with 



(85) 



(86) 



Applying here relation (|2lj) results in T/T = => T = 
const. The force projection which gives rise to such a 
behavior is 



u i F i m-(E-m)H 
A different special case is H 

^i9ik = -^-n-it ■ 



mUiF 1 



E 2 



m 2 ) H 



(87) 



equivalent to A = or T — 3H . According to Eq. (|2£ 
a production rate T = 3H is equivalent to an effective 
viscous pressure n = — (p + p). Via Eq. (|J) this implies 



8 



a constant energy density. Consequently, the case T 
3H corresponds to 



H 

H 2 







H = const 



a oc exp [Ht] , (88) 



i.e. exponential inflation. Both the CKV case on which 
the dynamics ( |S4|) is based and the case (^6|) correspond 
to constant ratios of F/3H, i.e., to constant values of A. 
Constant values of A (together with w = const) give rise 
to a fixed expansion dynamics of the scale factor, (either 
a specific power law or exponential inflation for the spe- 
cific choices above) which does not change in time. Since 
the early days of inflationary cosmology the graceful exit 
problem, i.e., the question of how to terminate a phase 
of accelerated expansion has played an essential role. A 
mechanism is necessary which induces a transition from 
accelerated to decelerated expansion. In the "standard" 
inflationary picture this mechanism is "reheating" , a non- 
equilibrium period during which the inflaton field decays 
into "ordinary" matter through a complicated multi-level 
process with copious production of particles (see p9"H32t ) . 
Obviously, a setting in which A is assumed constant, does 
not provide us with a transition between accelerated and 
decelerated expansion. On this basis it is impossible to 
obtain a cosmological scenario with a changing effective 
equation of state for the cosmic medium. A constant A 
may be considered realistic at most piecewise. To study 
more general situations in wich A is allowed to depend 
on time we proceed as follows 33 . We replace H in Eq. 

dH 

H = —a = H Ha . 
da 

where H' = AH /da. The resulting equation for the Hub- 
ble rate is 



H' _ 3 1 + w 
AH ~ ~2 a 



(89) 



For a constant A corresponding to T oc H it integrates 
to 



H _ /a \l( 1 + lu ) A 
Wo ~ \~a) 



(A = const) , (90) 



and 



a oc *3(i+«oa 



(A = const) 



(91) 



For A = 1/2 and w = we recover, of course, the be- 
havior (|84|). The case A = corresponds to exponen- 
tial expansion. If the linear dependence of the particle 
production rate on the Hubble parameter is replaced by 
a quadratic one, i.e., T oc H 2 , the quantity A will no 
longer be a constant. A dependence of the type T oc H 2 
was investigated in the context of matter creation due to 
cosmological vacuum decay [|4|3^]. It is equivalent to a 
particle production which is proportional to the energy 
density. With the ansatz 



r h 

3H ^TL 



(92) 



where the quantity H e = H (a e ) is the Hubble parame- 
ter at some fixed epoch with a — a e , integration of the 
differential equation j89|) for w — 1/3 yields 



H 



:H P 



(93) 



H changes continously from H = 2H e , equivalent to a oc 



exp [Ht] at a = 0, to H 



2 , equivalent to a oc t 1 / 2 , 



the standard radiation-dominated universe, for a ^> a e . 
For a < a e we have a > 0, while a < holds for a > a e . 
The value a e denotes that transition from accelerated 
to decelerated expansion, corresponding to H e = —H 2 . 
While a linear dependence of T on H corresponds to a 
power-law solution for the scale factor, we have now f3q] 



t = t P 



1 



4H e . 



In 



(94) 



Our point here is to obtain a microscopic realization of 
such a scenario in terms of a self-interacting one-particle 
force of the type @. From Eqs. @ and @ we find 
that for p = nT, p = 3nT, which is the ultrarelativis- 
tic limit of the equation of state (|50|), the linear part 
UiF*p a of the force projection vanishes. (Since the ultra- 
relativistic limit corresponds to m — * we use here and 
in the following a refedinition of the force according to 
mF l -> F l ). With the Hubble rate @ we obtain that 
the quantity A is given by 



A 



(95) 



For a — * we have A — > 0, i.e., T — * 3H, while for a ^> a e 
we find A — ► 1, i.e., r — * 0. This allows us to obtain the 
corresponding quadratic force term. The result is 



-EH 



(96) 



An internal interaction, described by this component 
of an effective one-particle force F l which is self- 
consistently exerted on the microscopic constituents of 
the cosmic medium, realizes the deflationary dynamics 
M). Eqs. © and @ specify to 



2H 



9ik 



a z + at 



:U t U k 



(97) 



where H is given by Eq. (|93|). The limit a — > corre- 
sponds to 



f 2H u 

£ u a gik — fHk 

T 1 



(a< a e ) , (98) 
while for a > a e we reproduces the CKV condition 
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p 2H 

<£ u„ Qik — —^r 



T 



Qik 



(a > o e ) 



(99) 



Formul a (|97[) continuosly mediates between an initial be- 
havior (|9q) and a final conformal symmetry (99), equiv- 
alent to a change between an initial phase with exponen- 
tial expansion where A = and a final stage in which 
the universe is characterized by a standard global equi- 
librium with A = 1, i.e. T = 0, for relativistic matter. 
Asymptotically, the vector u a /T is a CKV of the metric 
cjik- This brings us back to the question, to what ex- 
tent a modification of the conformal symmetry (e.g. in 
Eq. ((9^)) represents a symmetry again. To solve this 
problem we introduce the "optical metric" (cf. [p^|-|l^|) 



9ik = 9ik + 



nt 



UiUk 



1 



a e 

-§ = A ~ 



(100) 



where n r plays the role of a refraction index of the 
medium which in the present case changes monotonically 
from 7i r — ► 00 for a — > to n r — > 1 for a»8 e . For a = a e 
we have n re = \[2. Optical metrics are known to be help- 
ful in simplifying the equations of light propagation in 
isotropic refractive media. With respect to optical met- 
rics light propagates as in vacuum. Here we demonstrate 
that such type of metrics, in particular their symmetry 
properties, are also useful in relativistic gas dynamics. 
Namely, it is easy to realize that Eq. ( |97|) is equivalent 
to 



£ «« 9ik 



2H 



(101) 



The quantity u a /T is a CKV of the optical metric gik- 
This clarifies in which sense the modification of the con- 
formal symmetry is again a symmetry. According to Eq. 
( |100| ) the introduction of an optical metric corresponds to 
a mapping of the non-equilibrium contributions (A^l) 
onto a refraction index n r ^ 1 of the fluid. This map- 
ping is such that it associates the same conformal sym- 
metry to cjik which characterizes gik in the equilibrium 
case A = 1 with n r = 1 and gik — gik- The property of 
u a /T being a CKV during the entire process allows us 
to precise the nature of the fluid dynamics which realizes 
the deflationary scenario (|93|): The transition from a de 
Sitter phase to a FLRW period may be regarded as a spe- 
cific non-equilibrium configuration which microscopically 
is characterized by an equilibrium distribution function 
and macroscopically by the conformal symmetry of an 
optical metric with a time dependent refraction index of 
the cosmic medium. 



V. DECAYING VACUUM MODEL 

Combining the Friedmann equation ( |si"|) with the Hub- 
ble parameter (93h, we obtain the energy density 



ml 

2vr 



(102) 



where we have replaced G by the Planck mass mp ac- 
cording to G = 1/mp. Initially, i.e., for a <C a e , cor- 
responding to t — > —00, the energy density is constant, 
while for s»a e the familiar behaviour p oc a~ 4 is recov- 
ered. The temperature behaves as T oc p 1 / 4 . The same 
dependence is obtained for the particle energy E, such 
that E/T — const, consistent with Eq. (p0|) . The evolu- 
tion of the universe starts in a quasistationary state with 
finite initial values of temperature and energy density 
at a = and approaches the standard radiation domi- 
nated universe for a 3> a e . This scenario may alterna- 
tively be interpreted in the contex t of a decaying vacuum 
(cf. |3jj3^]). The energy density ( |102| ) may be split into 
P = P ( „) + P( T) where 



P(r) 



ml 
— -i 

ml 









\ al 1 




a 2 + a 2 e 



(103) 



The part p.. is finite for a — > and decays as a 6 for 
a a e , while the part p (r) describes relativistic matter 
with p.. — > for a — ► and p {r) oc a~ 4 for a ^> a e . The 
energy balances are (A = v, r) 



[p w 



r (A) [/>(, 



with 

r 



m V 2 a 2 



m 



2 a? 



The equation for p.. may be written as 

p w +3H(p M +p w +n_) =0 



where 

n, 



m 
1 



(104) 



(105) 



(106) 



1^ 
2 a 2 



This corresponds to an effective equation of state 



P. 



n, 



(107) 



(108) 



Although we have always = p M /3, the effective equa- 
tion of state for a — > approaches P, } = —p, v) ■ Effec- 
tively, this component behaves as a vacuum contribution. 
For a S> a e it represents stiff matter with P ( . = p.. . The 
radiation component may be regarded as emerging from 
the decay of the initial vacuum according to 



(109) 



10 



One may introduce a radiation temperature T. , charac- 
terized by a dependence T (r) oc P^ 4 - Similar thermo- 
dynamical consideration as those in sections II and III 
for the temperature T are possible for T, . . In particu- 
lar, one may consider the Lie derivative of the metric gik 
with respect to = u a /T {r) . Restricting ourselves here 

to the isotropic case, we have, analogously to Eq. (|20|), 



where 



T 

A - t±± 

(r) 



ik 



(110) 



(111) 



Along similar lines as in the one-component case we ob- 
tain [cf. Eq. @] 



£ 



Ti 



Ti 



(r) 



2H 



Ti 



(r) 



9ik 



2 a 2 



;UiU k 



(112) 



The quantity £° is a CKV of the optical metric g\ k , 



£ 



f (r) 9ik 



(r) 



2H 

T7, 



_(r) 



(113) 



a simple gas which is possible only for ultrarelativistic 
matter is no longer a singular case but a well defined 
interaction- free limit of a more general concept, (ii) It 
admits a simultaneos treatment of equilibrium states and 
a specific class of non-equilibrium states. The latter are 
characterized by an equilibrium distribution function of 
the gas particles (generalized equilibrium) although they 
may be far from equilibrium. 

We have established a microscopic realization of a de- 
flationary scenario of the early universe according to 
which the cosmological evolution started with a de Sit- 
ter period, followed by a smooth transition to a stan- 
dard radiation dominated FLRW phase. This deflation- 
ary dynamics represents a specific, exactly solvable non- 
equilibrium configuration of a self- interacting gas. Map- 
ping the non-equilibrium contributions onto an effective 
refraction index of the cosmic matter, the deflationary 
transition appears as the manifestation of a timelike con- 
formal symmetry of an optical metric in which the refrac- 
tion index changes smoothly from a very large value in 
the de Sitter period to unity in the FRLW phase. 
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where 



_(r) 

9 lk = 9ik 



1 - 



UiUk 



with 



,2 1 



1 + 



3/2 



(114) 



(115) 



This demonstrates that the dynamics of the radiation 
com pone nt is governed by a separate conformal symme- 
try ( |113| ) which, h owe ver, is closely related to the con- 
formal symmetry ( |l0l| ) for the med ium as a w hole. The 
corresponding optical metrics (114) and ( |100| ) have dif- 
feren t eff ective refractive indices which are connected by 
Eq. (pl!|). 



VI. CONCLUSION 

We have presented a theory of self-interacting gaseous 
fluids which systematizes and generalizes earlier in- 
vestigations concerning the relation between (non-) 
equilibrium and symmetry properties of the cosmic sub- 
stratum and the dynamics of the early universe. The 
general features of this approach are: (i) It provides an 
extension of the global equilibrium concept of relativis- 
tic gas dynamics. The standard global equilibrium for 
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